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A dynamic model for the simulation of gas micro combined heat and power devices (iulCHP boilers) has 
been developed in order to assess their energy performances. From a literature review and experimental 
investigations, the new model is designed with the aim of limiting the number of parameters which need 
to be easily accessible in order to be suitable with annual building energy simulations. At first, this paper 
focuses on the description of the iulCHP boiler which has been tested, on the development of the test bench 
and on the experimental results. Then, it focuses on the model description, on its parameterization and 
on its validation. The modelling approach is based on an energy balance on the device and on empirical 
expressions for the main inputs and outputs. The model characterizes the (jlCHP boiler behaviour for 
different inlet water flow rates and temperatures. The dynamic phases with the start-up and cooling 
phases are taken into account. Finally, the models for the Stirling engine and the additional boiler are 
limited respectively to 28 and 24. Further experimental investigations led to simplify the iulCHP model 
from 28 to 17 parameters without reducing the accuracy. 

© 2014 Elsevier B.V. All rights reserved. 


1. Introduction 

(jlCHP devices are high efficiency solutions for heat and electric¬ 
ity production which can reduce their environmental footprint by 
decreasing C0 2 emissions and by converting more primary energy 
than separate electricity and heating productions 1-4 . Moreover, 
piCHP devices can be a solution for demand-side management to 
relieve the power grid and to avoid using high carbon power gen¬ 
eration during the electrical demand peaks as shown by Vuillecard 
et al. 5]. CHP is already widespread in the industry and the district 
heating networks with a power range of MW. piCHP corresponds to 
low power CHP (less than 15 l<W el ) in order to be suitable for resi¬ 
dential or commercial heating and power needs. Because of the low 
power output of puCHP devices, the technologies which are used are 
different from centralized power plants [6 . The electrical efficien¬ 
cies are also lower but the global efficiencies are slightly higher. 
The main conversion processes which are used for piCHP applica¬ 
tions are: micro turbines, Stirling engines, Ericsson engines, steam 
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engines (Rankine cycle), internal combustion engines and fuel cells 
[7]. Among these technologies, the study focuses on the gas Stir¬ 
ling ptCHP boilers which are increasing on the market (Table 1). 
piCHP boilers consist in a Stirling engine sized to cover the heat and 
power base loads of an accommodation and an additional gas con¬ 
densing boiler to cover the peak heating demands (space heating 
and/or domestic hot water (DHW)). The excess electricity produc¬ 
tion is exported on the public grid. However, their introduction is 
slowed down due to their high investment costs in countries where 
electricity prices are low. The development of models suitable 
for use in building energy simulations is required for the assess¬ 
ment of these piCHP systems, notably in life cycle costs. Indeed, 
the potential growth of this technology on the market is highly 
dependent on investment costs, feed-in tariffs and investment sub¬ 
sidies. 

The piCHP boilers are generally controlled by the heating 
demand in order to maintain a high overall efficiency and to ensure 
the thermal comfort as a priority. The optimization of piCHP solu¬ 
tions including the sizing of the auxiliary components such as 
the thermal storage requires the development of models able to 
represent correctly their performance, especially during the start¬ 
up/cooling phases. 
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Nomenclature 

c 

heat capacity, J K -1 kg -1 

^ air 

air excess, - 

GHV 

gross heating value, J kg -1 

H 

enthalpy, JmoU 1 

H 

enthalpy flux, W 

LHV 

low heating value, J kg -1 

Lv 

latent heat, J kg -1 

M 

molar mass, kg mol -1 

777 

massflow rate, kgs -1 

[MC] 

thermal mass.JK -1 

P 

pressure, Pa 

P 

power, W 

Q 

heat flux, W 

t 

time, s 

T 

temperature, °C 

lUA ] 

exchange coefficient, WK -1 

^air 

combustion air volume, N m 3 N m -3 

V m 

molar volume, N m 3 mol -1 

V 

volume flow, N m 3 s -1 

Greek symbols 

AT 

temperature difference, K 

At 

time step, s 

X 

molar fraction, - 

T] 

efficiency, % 

P 

density, kg N m -3 

X 

time constant, s 

Subscripts and exponents 

0 

reference 

air 

combustion air 

amb 

ambient 

aux 

auxiliary 

cond 

condensation 

cool 

cooling phase 

cw 

cooling water 

el 

electrical 

electro 

electronic 

end 

end 

exh 

exhaust gas 

fuel 

fuel 

gross 

gross part 

HX 

heat exchanger 

i 

inlet 

int 

internal 

lat 

latent 

loss 

skin heat losses 

max 

maximum 

net 

net part 

nom 

nominal 

0 

outlet 

P 

power 

PV 

process value 

Q 

heat 

sen 

sensible 

SP 

set point 

start 

start-up phase 

stop 

stop phase 

th 

thermal 

tot 

total 


Besides, the model of |xCHP boiler to be developed for building 
energy simulation tools should meet the following requirements: 

• Short computation time in accordance with annual building 
energy simulations 

• Accuracy in steady states and transient (start-up and cooling 
phases) conditions 

• Easiness to parameterize 


2. Modelling approaches 

Many Stirling engine models have been developed for the 
description of the engine internal phenomena, and not for CHP 
applications coupled with the building. The usual time steps for 
accurate thermodynamic models (from the order of 10 -3 to 10 -6 s) 
are very different in comparison to building simulations (from 
15 min to 1 h). Moreover, to correctly model an engine, many 
features have to be known like heat exchanger surfaces, work 
pressures, geometry, etc., which are generally unknown. Many 
thermodynamic models on Stirling engines exist: Schwendig et al. 

[8] modelled the global Stirling cycle, Kongtragool and Wongwises 

[9] introduced the dead volumes effects. Organ 10] first pro¬ 
grammed CFD models of Stirling engines [11. Urieli and Berchowitz 
[12] worked on a simplified thermodynamic model in order to find a 
compromise between accuracy and computational time. However, 
this last model still requires many data, mainly on the regenerator 
and piston geometry. As a consequence, these detailed models do 
not match to the problem of coupling with the buildings. 

Beausoleil-Morrison and Kelly 13 have undertaken a literature 
review on models of mini-and p>CHP tools for dynamic simulation 
of buildings. Among the models identified in the literature review, 
Beausoleil-Morrison and Kelly mention McRorie and Underwood 
[14] whose model would suppose substantial effort to be integrated 
in building simulation codes, Kelly 15 who models a very specific 
two cylinders diesel engine, and the model of Pearce et al. [16] 
which does not take into account seasonal variations of efficiency. 
The authors 13] conclude that none of these approaches is suitable 
for coupled modelling to the building, since these models cannot be 
used to study control strategies and thermal storage which play an 
important role in the integration of p>CHP technologies in buildings 
[17]. 

Furthermore, “black box” models have been developed. Thiers 
et al. 18] have fitted a quasi-static model of a wood pellet Stirling 
engine from a measurement campaign conducted on the Sun- 
machine. The start-up and cooling phases are well characterized. 
This model is strictly limited to the tested product and to the 
achieved test conditions (source temperatures, water flow rates, 
fuel characteristics, heating load). Another black box generic model 
is proposed in the French thermal regulation which uses a constant 
annual efficiency [19] assessed by a short time experimental test. 
The efficiency is independent on any operating conditions. Conroy 
et al. [20] made a numerical model of the gas Stirling engine p^CHP 
Whispergen based on field tests. The accuracy of the measurements 
is not given and the fitting only depends on time without physical 
meaning. Roselli et al. [21 also studied the energy performances 
of 8 |jiCHP devices including 3 power modulating devices. Steady 
states have been characterized but no numerical model has been 
implemented. 

The model of Annex 42 of the International Energy Agency [13] 
is the closest to our set objectives. It is adapted for the implemen¬ 
tation of p>CHP models (internal combustion engine and Stirling 
engine) in dynamic simulation piece of software. This model uses 
2 thermal capacitance nodes. An optimization procedure is used 
to identify the 91 parameters but the uniqueness of the solution 
(identified parameters) is not guaranteed. Lombardi et al. [22] have 
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Fig. 1 . Test bench (left) and |jiCHP (right).* 1: exhaust gas outlet, 2: condensation heat exchanger, 3a & b: burners, 4a: head of the Stirling engine, 4b: Stirling engine, 5: fan, 
6a & b: gas suppliers, 8: water pump, 9: DHW heat exchanger). 


simplified this model by implementing a linear difference model 
based on linear coefficients and reference parameters which are 
determined by a few experimental tests. Borg and Kelly [23] also 
developed a dynamic model for tri generation studies, based on 
inflows and outflows of the system. They also underline the rel¬ 
evance to take into account the dynamic behaviour of energetic 
systems to obtain realistic results in comparison with static models. 

This paper presents a generic model of |jiCHP boilers, based on 
experimental investigations on a gas Stirling engine p,CHP system, 
with a few parameters easily reachable by simple experiments or 
published characteristics. Each parameter is determined indepen¬ 
dently without global minimization. 


3. Experimental investigations 

3.2. Experimental set up 

The tested unit is a gas Stirling engine p^CHP from De 
Dietrich/Remeha called Hybris Power including an additional con¬ 
densing gas boiler. Hereafter, both the Stirling engine and the 
additional boilers are separately studied. The test bench is designed 


to measure the main values involved in the energy balance. The 
bench gets several operation modes: 

• Instantaneous heating (heating and DHW) and electric produc¬ 
tion. 

• Semi stored heating: a part of the heat production is stored in a 
storage tank of 1301. 

Figs. 1 and 2 present the bench and the metrological means: 
The metrological means (Table 2) have been selected to provide 
acceptable uncertainties for the fuel consumption, electrical and 
thermal outputs, which are comparable with other experimental 
trials [24-26]. 

The thermal power is calculated by the thermal energy transport 
equation (Eq. (1)): 

QhX = P(Tcw,i ) ' Vcw • i c cw{Tcw,o ) • Tcw,o ~ c cw(T CW j) • T cw j ) (1) 

The thermal properties of the water are obtained from standard 
tables. The fuel power is calculated by measuring the gas volume 
flow converted to normal conditions of pressure and temperature 
and by using the GHV of the gas daily given by the supplier. 

Pfuel = m/ueiGHV ( 2 ) 



Fig. 2. Hydraulic principle scheme. 
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Table 1 

Existing gas |jiCHP prototypes. 


Manufacturer 

Apparatus 

Power 


LHV efficiency 




Electric, kW e i 

Thermal, kW t h 

Electric, % 

Thermal, % 

Gas Stirling engine (gas SE) 

Enatec 

Infinia 

1 

6.4 

12.5 

80 

SOLO 

Stirling 161 

2-9.5 

8-26 

22-24.5 

65-75 

Disenco 

Inspirit 

3 

15 

16 

76 

Baxi 

Ecogen 

0.3-1 

3.7-25.3 

16 

83 

Viessmann 

Witowin 300-W 

1 

3.6-20 

15 

82 

WhisperGen 

EU1 

1 

7.5-14.5 

11 

84 

De Dietrich Remeha 

Hybris Power 

1 

3-23.7 

17 

85 

Senertech 

Stirling SE 

1 

3-23.8 

14 

77 

Sunmachine 

Gas 

1.5-3 

8-15 

25 

65 


: available. 

: introducing phase. 


The electrical power is measured by an energy metre which 
measures the incoming electricity consumed by the auxiliaries (fan, 
pumps, electronic card) and the out coming electricity produced by 
the Stirling engine. Uncertainties are calculated according to the 
propagation law of uncertainties [27,28]. 

3.2. Experimental results 

32 A. Stirling engine in steady state 

The main results have been plotted in Figs. 3-7. The tests show 
that the inlet water temperature has a significant influence on the 
electrical production, the condensation gains, the fuel power and 


the exhaust gas temperature. The variations of the thermal power 
are low and are located in the range of the uncertainties (Fig. 3). 

The theoretical Stirling cycle efficiency depends on the temper¬ 
atures of the cold and hot sources. The inlet water temperature 
represents the cold source and the burner supplies the hot source 
of the engine. Low inlet water temperatures increase the temper¬ 
ature gap between the two sources and increase the theoretical 
efficiency. For the condensation gains, the heat exchanger surface 
has to be the coldest as possible to condensate the steam from 
exhaust gas. Lowest inlet water temperatures in the heat exchanger 
can achieve significant heat gains. The mass flow rate also has an 
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Fig. 3. Thermal power of the Stirling engine. 


Fig. 5. Gross electrical power of the Stirling engine. 



25 30 35 40 45 50 55 60 65 T cw4 [°C] 



Fig. 4. Fuel power of the Stirling engine. 


Fig. 6. Condensation gains of the Stirling engine. 
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Metrological means. 
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Type 

Notation 

Metrologic means 

Error 


Consumed fuel 

V fuel 

Pulse counter 

±1.5% 


Fuel pressure 

Pfuel 

Hand pressure metre 

±1 mbar 

Fuel 

GHV 

GHV 

Gas supplier data 

±0.1 kWh N m -3 


Temperature 

Tfuel 

PtlOO probe 

±2 K 


Fuel power 

Pfuel 

Calculated 

±2.5% 


Inlet water temperature 

Tcw.i 

Pt500 probe 

±0.2 K + 0.05% 

Heating output 

Outlet water temperature 

Volume flow 

T cw ,o 

m cw 

Pt500 probe 

Pulse counter 

±0.2 K +0.05% 

±0.15 Lmn - 1 +2.5% 


Thermal energy 

Qhx 

Calculated 

±4.5% 

Electric output 

Electric power 

Pnet 

Watt metre 

± 1 % 

Exhaust gas 

Temperature 

Gas analysis: O 2 , NO x , CO, CO 2 , THC or CH 4 

Pexh 

PtlOO probe 

Combustion analyser 

±0.3 °C 


influence because it plays on the mean temperature in the heat 
exchanger. The greater the flow, the lower the temperature differ¬ 
ence between the inlet and the outlet is large at constant power. 
Low mean temperatures in heat exchangers lead to higher perform¬ 
ances on the electrical efficiency and on the condensation gains. 

These trends are confirmed by experimental investigations. The 
thermal output can be almost considered as constant, there is 
no modulation here. To maintain the thermal output as constant 
(Fig. 3), the combustion power increases with the inlet water tem¬ 
perature to compensate the latent heat losses which are greater 
(Fig. 4). The electrical power increases with low inlet water temper¬ 
atures and high water volume flow rates (Fig. 5). The condensation 
gains increase with low inlet water temperatures and high volume 
flow rates (Fig. 6). The exhaust gas temperatures increase with high 
inlet water temperatures and slightly decrease with higher water 
volume flow rates (Fig. 7). 

For the electrical power, the condensation gains and the exhaust 
gas temperature, the graphs use the volume mass flow for the x- 
axis to highlight the volume flow rate impact on these quantities 
which is less significant than the inlet water temperature. 

The condensation gains are assessed at the outlet of the device 
by measuring manually the condensates mass during a given time 
that gives a mean condensate mass flow rate. The sensible heat is 
not taken into account, only the latent heat is used to calculate the 
power recovery made on the condensation of the steam in exhaust 
gas. 

3.2.2. Additional gas boiler in steady state 

The same tests have been carried out on the additional boiler. 
The additional gas condensation boiler can modulate its thermal 
power. Figs. 8-11 give the energy performances of the additional 
boiler. The tests use a water flow rate range different from the 
Stirling engine tests because the involved powers are higher. To 
avoid high temperature differences and too high outlet water 
temperatures, the nominal volume flow rate of the Stirling engine 



(which also is the maximum value of the tested range for the 
Stirling engine) is used as the minimal volume flow rate for the 
additional boiler. The boiler thermal and fuel outputs are more 
sensitive to the inlet water temperature than the Stirling engine 
(Figs. 8 and 9). The volume flow rate also has a significant impact, 
linked to the internal control of the device. These two quantities 
decrease with high inlet water temperatures and low water 
volume flow rates. The additional boiler follows the same trends 
than the Stirling engine for the condensation gains and the exhaust 
gas temperatures (Figs. 10 and 11). 

3.2.3. Transient states 

The different transient phases have been studied. The results 
have been normalized in order to identify the general behaviour 
for any configurations and power levels (Figs. 12-17). The experi¬ 
mental data are fitted by simple laws like exponential or linear laws. 
Different tests have been carried out by changing the inlet water 
temperature and the volume flow rates. The thermal behaviour 
for the Stirling engine is more variable according to the test 



25 30 35 40 45 50 55 60 65 T cwi [°C] 

Fig. 8. Thermal power of the additional boiler. 



25 30 35 40 45 50 55 60 65 T cwi [°C] 


Fig. 7. Exhaust gas temperature of the Stirling engine. 


Fig. 9. Fuel power of the additional boiler. 
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Fig. 13. Normalized thermal power during the cooling phase of the Stirling engine. 


Fig. 10. Condensation gains of the additional boiler. 



12 14.5 17 19.5 22 V cw [Lmtl' 1 ] 

Fig. 11. Exhaust gas temperature of the additional boiler. 


conditions. The internal control and the start-up gas over-supplies 
make the thermal power variable during this phase. However, to 
keep a consistency of the model, the general trend of the thermal 
output can be fitted by an exponential law in this phase by main¬ 
taining an energy equivalency (Fig. 12). The cooling phase is more 
stable with a similar exponential drop trend for different test con¬ 
ditions (Fig. 13). The electrical production has a similar behaviour 
during the start-up and cooling phases for any inlet water tem¬ 
peratures or water volume flow rates (Figs. 14 and 15). The drop 
is linear for the electrical production during the cooling phase 
(Fig. 15). The behaviour of the thermal power during the start-up 
and cooling phases for the additional boiler is stable for each test 
condition and exponential laws are used to fit the experimental 



Fig. 14. Normalized electrical power during the start-up phase of the Stirling engine. 



Fig. 15. Normalized electrical power during the cooling phase of the Stirling engine. 

data (Figs. 16 and 17). These fittings are used later in the model 
development. 

The tested p>CHP has only one electrical wire for outgoing and 
incoming electricity. It is impossible to dissociate the gross and net 



Fig. 12. Normalized thermal power during the start-up phase of the Stirling engine. 
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0 300 600 900 1200 1500 1800 Time [s] 

Fig. 16. Normalized thermal power during the start-up phase of the additional 
boiler. 

production in steady state. So, the cooling phase is used to esti¬ 
mate the net electrical power: after the stop of the Stirling engine 
electrical production, only the pump, the fan and the electronic 
consumption are operating. Three different cooling phases are pre¬ 
sented in Fig. 13. These tests also lead to the determination of the 
cooling phase duration. 

Up to At=90s from the beginning of the cooling phase, the 
engine keeps producing electricity with a linear drop (Fig. 15). A 
desynchronisation appears when the head temperature decreases 
below 170°C (without electrical production, the net electrical 
power becomes negative). Up to At=750s, the fan and the pump 
keep working to ensure the engine cooling until a head tempera¬ 
ture of 130 °C. The measured consumed power is between 40 and 
50 W (the resolution display is 10 W). A hand-held wattmeter has 
been used in order to get more accurate data of 44 W for the global 
auxiliary consumption. From At = 750 s, the p,CHP is stopped; only 
the electronic power is consumed for controllers. Manual measure¬ 
ments gave an electronic power of 4 W. 

4. Dynamic model development and identification 
procedure 

The dynamic model and the identification procedure are based 
on four operation periods: start-up phase, steady state, cooling 
phase, and stop phase. Three kinds of parameters are distinguished: 
steady state parameters, transient parameters and operation logic 
parameters which depend on the own control strategy of the unit 
and which are integrated in the transient parameters (Appendix A). 

The requirements of the developed model to be implemented 
in building energy simulation tools are summarized: 

• The time step has to be equivalent to building simulation time 
steps: some minutes. This point excludes accurate thermody¬ 
namic models that need tiny time steps (some seconds or less). 



Fig. 17. Normalized thermal power during the stop phase of the additional boiler. 


• The number of parameters to identify should be limited (an easy 
identification based on experiments is aimed where each param¬ 
eter keeps a physical sense). 

• The model be adaptable to other p>CHP units. 

• The model should be representative of the p,CHP units perform¬ 
ances (in particular, thermal and electrical transient behaviours 
of |xCHP during the start-up and cooling phases). 

In the following, the modelling is split in two parts according to 
the operating modes: 

• Steady state. 

• Transient states (start-up, cooling, stop). 

4.1. Steady state 

In steady state conditions, the electrical power, the fuel con¬ 
sumption and the heating power are assessed by correlations 
dependent on the inlet water temperature T cw , f , the water flow rate 
rricw , and the thermal load in the case of modulating p>CHP as pro¬ 
posed by Lombardi et al. [22] and Andlauer [29 . C0 2 emissions 
are deduced from the fuel consumption. The steady state perfor¬ 
mance is characterized by linear functions for the thermal power 
and the electrical power outputs and the fuel thermal input. The 
linear coefficients are identified from experimental tests in steady 
state conditions. 

In the following, the equations are written for non-modulating 
engines (ON-OFF operation) since the tested p^CHP cannot modu¬ 
late his load. In case of engines able to operate at variable capacity, 
the correlations of the net electric power and the fuel consumption 
should integrate a load factor on the thermal power accounting for 
performance degradation/improvement with the modulation. 

The expressions of the delivered thermal output, fuel input and 
gross electrical power, are defined in functions of the nominal con¬ 
ditions of flow rate and temperature. As shown in Figs. 3-5, and as 
proposed by Lombardi et al. [22], the equations can be written as a 
linear combination of water flow rate and inlet water temperature: 

P fuel = PfueT + (W - C|?) + ^ (mew - <T) (3) 

Qhx = Qj\°x m + c ( Tcw,i - TZ m i ) + d (mew - (4) 

Pgross = P™™ + * (W " T Z”D ^ ~ (5) 

If reference characteristics for the temperature or mass flow rate 
are given by the manufacturer, it is suitable to use it as nominal 
conditions. Otherwise, a method is given here to determine rele¬ 
vant nominal conditions. Concerning the nominal electrical output, 
the manufacturer value is used. For the inlet cooling water temper¬ 
ature, it is relevant to use a low temperature in order to optimize 
the condensation gains: 30 °C for example. The nominal mass flow 
rate has to match the temperature difference given by the man¬ 
ufacturer when the |jiCHP operates without the additional boiler. 
Here a temperature difference of 5 K has been fixed that led to a 
nominal mass flow rate of 16 kg mn _1 . 

In order to parameterize Eqs. (3)-(5) according to the cooling 
water mass flow rate and temperature, several tests should be 
achieved. The experiment plan consists in changing one parameter 
on the full test range while keeping the others in nominal condi¬ 
tions. The proposed ranges for the Stirling engine are compiled in 
Table 3. 

Net electric power corresponds to the usable part of electricity 
produced by the engine. Another part of the production is auto con¬ 
sumed by the auxiliaries (pumps, fans) P aU x and the electronic card 
Pelectro We define: 

Pgross = Pnet + Paux + Pelectro 


(6) 
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Table 3 

Tests list for the parameterization procedure and validation for the Stirling engine. 


rhcw 


0.4 rh™ m 


0.6 rhZ m 


0.8 mZ m 


j^nom 

,n cw 


Tcw,i 

j'nom 

cw,i 

T n °m + 1Q 
cw,i 

T nom + 20 

CW,l 

T nom + 30 

CW,l 


Electrical consumptions linked to the operation of the auxiliary 
devices like the fans and the pumps can be determined by differ¬ 
ent ways: either the p>CHP is equipped by two different electrical 
wires (one for the incoming electricity and one for the outgoing 
electricity) or these consumptions can be deduced during the cool¬ 
ing phase where no electricity is produced but the auxiliary devices 
are still working. These values can also be found in manufacturer 
data. 

For each test, the mean value of the measurements during the 
steady state is taken into account. The steady state is reached when 
the following conditions are obtained: 


(T sp . - 1.5) < T pv . < (T sp . + 1.5) and 0.95 x P* p 

S CW I ' r\M i ' c~\m i / H/C 


CWJ 


CWJ 


gross 


< 


pPV 


gross 


< 


1.05 x P|P ss . 

As shown in Figs. 9 and 10, the same equations as for the Stirling 
engine can be used: 


p f _ pnom 

1 fuel aux ~ 1 fuel aux 


+ a' (T, 


nom 


CWJ 


T™"!) + b' (mcw-nv 


QHXaux = + c' (T, 


CWJ 


TZj) + d ' ( m C w-rhZ m ) 


(7) 

( 8 ) 


The parameterization procedure is similar to those used for 
IJiCHP. The identified parameters on the tested p,CHP boiler are 
given in Appendix A. 

Assuming the combustion is complete, the fuel mass flow rate 
can be assessed from the fuel thermal input as: 


Pfuel 

= GHV 



The model can be used for different types of fuel. For the model 
identification, the fuel mass flow rate should be measured and the 
gross heating value (GHV) of the fuel should be known or measured. 
The air excess should be also known. From the fuel mass flow rate, 
the fuel composition and the air excess, typical combustion calcu¬ 
lations are undertaken in order to assess the combustion reactants 
mass flow rates and the exhaust gas heat losses. In addition to 
the steady state characterization, the start-up and stop transient 
phases are to be modelled. 


4.2. Transient phases 



0 15 30 45 60 75 90 105 120 time [min] 

Fig. 18. Normalized electrical cycle efficiency according to the cycle duration. 


During the Stirling engine start-up phase, we distinguish three 
phases. At the beginning, the Stirling engine head has to be heated. 
In a natural gas engine, this heating phase is achieved by a gas 
burner with a fuel mass flow rate greater than nominal mass flow 
in order to quickly reach the suitable temperature. In the model, no 
overshooting is assumed and the fuel mass flow rate is equal to the 
nominal conditions. As soon as the head reaches the set tempera¬ 
ture, the free piston begins to move and the electrical production 
can begin. Synchronized with the network, the p,CHP produces the 
electricity. 

We assume a first order behaviour for the electrical power, heat¬ 
ing power and head temperature with phase shifts as experimental 
investigations show. 

• During Start-up 1, the head temperature rises with gas combus¬ 
tion in the burner. We assume heating production is insignificant. 

• During Start-up 2, the combustion continues and the heating pro¬ 
duction starts, except for the electrical production. 

• During Start-up 3, the combustion and the heating production 
continue and the electrical production starts. 


The heating and electrical productions during the start-up phase 
are assumed to increase as an exponential function using delay 
times and time constants (Eqs. (10) and (11)). 


o 


fort < Utan + AtS, 


start 


Qhx = < 


/ t - (t s tart + At^ arf ) \ 


nom 

!X 


1 - e 


start 


fort > t star t + A tf tan 


J 


Pnet — \ 


0 for t < tstart + At® art 

f “ ( tstart + Atf tort ) 

fort > t s tan + At s P Mr( 


( 10 ) 


( 11 ) 


P„T I 1 - e 


start 


The use of exponential curves for transient phases gives a good 
representation of the evolution during these phases even if the 
curves are more complex due to the desynchronisation of the power 
generator during the cooling phase and fuel overshooting during 
the preheating phase. 

4.2.2. Start-up transient 

The start-up phase differs according to the engine type. The fol¬ 
lowing approach can be only partly considered as generic since it 
has been developed for the specific case of a Stirling engine. For 
other technologies, the model has to be adapted. Indeed, a fuel 
overshooting for start-up can appear. Moreover, when the engine 
is not fully cooled down at the ambient temperature, the start-up 
is partial. Andlauer [29] shows that the partial start-up mode can 
be considered as a full start-up mode on the specific Stirling engine 
studied in this paper. 


These parameters can be identified by experimental investiga¬ 
tions (Section 3.2.3). Physical models exists by taking into account 
the inertia of the heat exchanger or the mass of water [30], but 
the goal here is to propose a simple approach of modelling which 
requires the less parameters as possible. The identified parameters 
are given in the Appendix A. The description of the modelling of the 
start-up phase is given in Appendix B. 

The modelling of the transient phase lets to take into account 
the temporal inertia of these systems which have to match the heat 
and power demands. The profitability of these systems is linked to 
the electrical efficiency but also to the duration of the cycle. Fig. 18 
gives the evolution of the ratio between the electrical cycle effi¬ 
ciency according to the cycle duration and the nominal efficiency 
in steady state (normalized electrical efficiency). It shows that to 
reach high-normalized electrical efficiencies (95% of the nominal 
value for example), the system has to operate at least 1 h. 
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4.2.2. Cooling transient 

When the combustion stops inside the p^CHP engine, the pumps 
and fans keep operating for a given period called “cooling”, in order 
to use the residual heat and to cool down the engine. Then, these 
auxiliary devices are switched off. The heating production drop is 
approximated by a first order transfer function, characterized by a 
time constant (Eq. (12)). A partial cooling mode could be taken into 
account but it is not implemented in the model developed. 




0 


Qhx 



t-t 


stop 


A 

'stop 


fort > tstop + At^op 
fort < t s top + At£op 


Pnet — 


-Paux Pelectro 

t — tstop 


fort > tstop + A 


stop 


pnom i _ 
1 gross 1 1 


A t P 

^ L stop 


Paux Pelectro fort < t s top + At5r 


stop 


( 12 ) 


(13) 


The time constants are to be identified on cooling phases 

(Figs. 13 and 15). 

The |jiCHP remains able to produce electricity thanks to the accu¬ 
mulated heat on the hot side of the engine. The electrical production 
lasts a short period compared to the pump and fan operation time 
for the engine cooling. The electrical output decreases until 0, but 
not necessary exponentially (Fig. 13). The proposed method con¬ 
sists in fitting the electric drop by a linear law (Eq. (13)). 

The cooling phase period, while the fan and the pump keep 
operating, depends on engine control: either the engine head tem¬ 
perature reaches a limit or a certain time is elapsed. The model 
uses a time limit, after which the cooling phase stops. This choice 
is adapted for the tested Stirling engine since the engine is non¬ 
modulating, so that the head temperature is always the same at 
the beginning of the cooling phase and thus the cooling time. 

During the engine stop phase, after the cooling phase, the fan 
and the pump are stopped. The engine temperature evolution is 
only governed by heat exchanges with the ambiance Q/ oss , the ther¬ 
mal, fuel and electrical powers are equal to 0. Only the electronic 
power P e i ec tro is consumed. 

The model parameters for the cooling phase of the tested p>CHP 
are presented in Appendix A. The modelling approach of the cooling 
phase is described in the Appendix B. 


4.3. Exhaust gas heat loss 


The model parameters for the condensation gains of the tested 
p,CHP are presented in Appendix A. 


4.3.3. Enthalpy flows 

The exhaust gas latent heat losses, that is to say the latent energy 
of the non-condensed fraction of steam are defined in Eq. (18). 

H% h = (m™ - m cond )L v (18) 

The maximum condensation mass flow rate is assessed as the 
mass flow rate of the steam in the exhaust gas based on the fuel 
composition and the fuel mass flow rate. The reactants and com¬ 
bustion products enthalpies are calculated fromJANAF tables [31] 
giving for each component a polynomial expression of the specific 
enthalpy. 

Hi = > prY, {[Hi{T) - Hi{To)]Xi] (19) 

i 

Finally, the exhaust flow rate becomes: 

P^exh = PI fuel T PI air ~ P^cond (20) 

With: 

m Qir = ihfuelVair (1 + e Qir )^ (21) 

v m 

4.4. Dynamic model 


If the main outputs of the model are given by the correlations 
presented above, the use of a dynamic model provides information 
on the behaviour of the p>CHP such as the repartition of heat losses 
and dynamic phases of the engine. The proposed model uses only 
one thermal capacitance for the whole p^CHP device. This assump¬ 
tion is made to facilitate the model parameterization; the start-up 
and cooling periods are each represented by only one time constant. 
A fictitious “global average internal temperature Tint is introduced 
which governs the heat losses between the p^CHP skin and the 
ambiance using a fictitious coefficient [DA] /oss . 

The energy balance can be written as: 


[MC] 


dTi 


int 


int 


dt 


— Pi fuel T H n j r -T Pfupj Pi 
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gross 


Qhx - Qi 


OSS 


risen 

n exh 


iflat 
n exh 


( 22 ) 


4.3.1. Exhaust gas temperature 

The exhaust gas temperature is modelled with linear corre¬ 
lations using the cooling water mass flow rate and temperature 
(Figs. 7 and 11), like the other main inputs or outputs. 


T exh =g(T cwJ -TZ m i ) + h (mew - mZ m ) + T; 


nom 

exh 


T _ (r f (T 'rnorrn , ur p^nom \ , ^nor 

1 exhaux — § \*cw,i ^cw,i^~^~^ l EUcw P^cw ) T j £ xh 


nom 
aux 


(14) 

(15) 


The model parameters for the exhaust gas temperature of the 
tested p,CHP are presented in Appendix A. 


Qloss = [ DA] i oss • (T int - T amb ) (23) 

Due to the conventional definition of T int , [DA] /oss is strictly rela¬ 
tive to its assessment. A simple formula giving T int based on thermal 
measurements (based on the inlet cooling water temperature and 
the hottest temperature measured on the hot side of the Stirling 
engine) is proposed by Andlauer and al 29,32 . The first order dif¬ 
ferential equation can be solved at each time step by using a finite 
difference method. 

The model is presented in a flow diagram in the Appendix C. 


4.3.2. Condensation gains 

The condensation gains are modelled by linear correlations. 
From Figs. 6 and 10, it can be defined as: 


rheond = max (0; m”™ + i (T Wii - T™) +j (mew - )) (16) 


m co ndaux = max (0;m 


pnom 


cond aux 


4 - i' (T ■ - juom\ 
^ 1 t 1 cw,i 1 cw, i ) 


+f (ihcw-mSD) 


(17) 


5. Simplified model 

A model which can extrapolate the system performances in 
different operation ranges even when the control strategy or the 
design changes is very useful and relevant provided that the num¬ 
ber of parameters is limited and their identification requires few 
experimental investigations [33 . The model of the Stirling engine 
needs to identify 28 parameters which can be estimated by experi¬ 
mental investigations as shown in Section 4. In this part, a reduction 
of the model parameter number is proposed in order to facilitate 
the model parameterization. 
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Table 4 

Comparison between the model and experiments during steady state for the Stirling engine. 



T °r 

29.8 

40.2 

50.1 

60.0 

29.9 

40.4 

50.2 

60.0 

29.8 40.5 

50.3 

59.9 


mew Is- 1 

8.5 

8.5 

8.5 

8.5 

13.2 

13.2 

13.2 

13.2 

16 16 

16 

16 


Pfuei 

-1.3 

-0.3 

-2.5 

3.2 

-0.1 

1.1 

-3.3 

0.3 

3.4 -0.5 

-2.4 

0.4 


CW 

-3.1 

0.1 

2.8 

-0.6 

-1.2 

1.2 

0.9 

-1.6 

2.7 -0.6 

-0.9 

-1.1 

Conrm ete method % 


Pgross 

0.6 

-0.9 

-1.1 

1.4 

1.1 

-0.8 

-1.2 

0.8 

0.7 -0.1 

-1.2 

0.6 


Texh 

-0.6 

-1.3 

-0.9 

3.4 

4.9 

1.6 

-2.4 

1.7 

1.3 3.4 

-0.9 

1.2 


Pfuei 

-2.7 

-1.7 

-3.9 

1.8 

0.4 

1.6 

-2.9 

0.8 

5.1 1.1 

-0.9 

2.0 


CW 

-2.9 

0.5 

3.5 

0.5 

-1.8 

0.9 

0.9 

-1.4 

1.7 -1.4 

-1.4 

-1.3 

SimDhhed method 


Pgross 

0.9 

-0.6 

-0.8 

1.7 

1.0 

-0.9 

-1.3 

0.7 

0.4 -0.5 

-1.5 

0.2 


Texh 

-8.3 

-5.6 

-3.2 

2.7 

0.4 

0.0 

-2.2 

3.3 

-0.9 3.6 

0.9 

4.0 


Table 5 

Comparison between model and experiments during steady state for the additional boiler. 



t . °r 

i CW,l> 

29.9 

40 

49.7 

59.1 

29.9 

40.2 

50.3 

61.1 

29.8 

39.9 

50.4 

57.4 


riw Is -1 

12 

12 

12 

12 

17 

17 

17 

17 

22 

22 

22 

22 


Pfuei aux 

5.9 

-3.8 

18.6 

16.3 

-1.5 

-4.9 

-9.3 

4.9 

7.6 

1.6 

-5.3 

-5.1 

Complete method (%) 

QhXoux 

0.3 

-8.8 

11.6 

5.4 

-1.1 

-1.4 

0.5 

2.6 

5.5 

1.0 

-2.4 

-7.7 


Texh aux 

2.5 

0.4 

-6.9 

2.8 

1.7 

-3.4 

-5.6 

3.4 

2.9 

-2.3 

-5.4 

5.9 


Table 6 

Comparison between model and experiments during transient states for the Stirling engine. 



T °C 

1 CW,l y ^ 

29.8 

40.2 

50.1 

60.0 

29.9 

40.4 

50.2 

60.0 

29.8 

40.5 

50.3 

59.9 


mew Is -1 

8.5 

8.5 

8.5 

8.5 

13.5 

13.5 

13.5 

13.5 

16 

16 

16 

16 

Start-up (%) 

Qhx 

12.0 

-10.3 

-9.2 

13.3 

-12.3 

-12.4 

8.7 

9.2 

-9.2 

10.0 

-7.3 

-6.3 

Pgross 

0.9 

-0.6 

1.1 

0.8 

-1.8 

1.1 

-2.0 

1.3 

-1.5 

0.4 

-1.0 

-0.3 

Cooling (%) 

■1 

1.2 

2.1 

3.2 

-1.2 

-1.8 

4.1 

-1.3 

2.0 

3.4 

1.1 

-2.7 

0.9 

Pgross 

5.3 

4.2 

-2.3 

3.7 

-1.0 

2.0 

-4.2 

1.7 

-6.1 

-2.1 

5.0 

2.4 


Table 7 

Comparison between model and experiments during transient states for the additional boiler. 



t • °r 

29.8 

40.2 

50.1 

60.0 

29.9 

40.4 

50.2 

60.0 

29.8 

40.5 

50.3 

59.9 


HTcw IS" 1 

12 

12 

12 

12 

17 

17 

17 

17 

22 

22 

22 

22 

Start-up (%) 

Qhx aux 

-8.2 

1.2 

2.7 

-5.0 

5.1 

-7.6 

3.0 

-5.8 

2.8 

-9.7 

-4.1 

5.6 

Cooling (%) 

Qhx aux 

6.6 

9.2 

-2.1 

2.2 

-11.9 

4.2 

-7.0 

4.6 

-6.0 

7.2 

-1.1 

7.0 


5.2. Reduction of the number of parameters 

Eqs. (3)-(5), (14), and (16) for the fuel power, the heating power, 
the electrical power, the exhaust gas temperature and the conden¬ 
sation gains can be simplified. As shown in Figs. 3, 4, 19 and 20, 
the influence of the cooling water mass flow rate can be neglected 
and the influence of the inlet cooling water temperature, which is 
more significant, can be characterized by a linear function. For the 
heating power, it can be simplified as a constant value (Fig. 3). 

The different equations become: 


Qhx = OUT (25) 

Pgross = + e (T CW J - C™) (26) 

T'exh = T cw j + AT (27) 

The last simplification consists to include all the heat losses 
(latent and sensible exhaust gas heat losses, skin heat losses, air and 
fuel enthalpy flows) in one heat loss source limiting the information 
provided by the model but avoiding an experimental characteriza¬ 
tion of the condensate mass flow rate. 


3 _ pnom i a (t 

fuel ~ 1 fuel + ^ 1 < 


CWJ 


Tnom\ 
1 cw, i' 


( ' 24 '* H/ oss = Pfuei - Pgross - Qhx 


(28) 




Fig. 19. Gross electrical power of the Stirling engine. 


Fig. 20. Exhaust gas temperature of the Stirling engine. 
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5.2. Use of manufacturer data 

With these simplifications, the model requires to identify only 
17 parameters. Some manufacturer data are generally available 
and could be used to identify several parameters. Main behaviour 
laws are linear so the modelling just needs two performance points 
to characterize the unit. It seems reasonable to get performance 
data from manufacturers for two operation points at two differ¬ 
ent temperatures, Ti and T 2 . The auxiliary and electronic powers 
of the pCHP and the air excess can be obtained from the manu¬ 
facturers. The remaining parameters such as time constants seem 
more difficult to obtain from manufacturers. Default values could 
be chosen or identified by an experimental test with a start-up 
and a cooling phase. The simplified model parameters are listed 
in Appendix A. The reduction of parameters is only applied on the 
Stirling engine. On the additional boiler, only the simplification on 
the exhaust gas and the aggregation of the heat losses are used, the 
mass flow rate has a strong impact on the outputs and cannot be 
neglected. 

6. Validation of the model and discussions 

The model with the two-parameterization methods (28 and 17 
parameters) is compared to experimental results in steady and 
transient states. The experimental measurements are carried out 
after more than 3h of steady state operation times. The tests 
are carried out for different inlet cooling water temperatures and 
cooling water mass flow rates. Table 4 gives the relative differ¬ 
ence between the two models of pCHP and the experimental 
results. 

The relative differences between the simplified and the com¬ 
plete parameter identification procedures show that the both 
methods are very suitable to correctly model the real behaviour 
of the unit. 

Table 5 shows the results of the comparison between the addi¬ 
tional gas boiler model and the experimental results. 

Tables 6 and 7 show the relative differences between the Stirling 
engine and the additional gas boiler models and the experimen¬ 
tal results for the transient states. The modelling method of the 
transient phases by time constants is very suitable to correctly 
model the real transient behaviour of the unit except for the ther¬ 
mal output during the start-up of the Stirling engine where the 
difference is higher at about 10% because of the more chaotic 
behaviour. 

The Stirling engine can work simultaneously with the additional 
boiler. In this configuration, the inlet water temperature of the addi¬ 
tional boiler matches to the outlet water temperature of the Stirling 
engine. The steady state and transient inputs and outputs are added 
by using the superposition principle. 

The efficiencies of this system decrease with higher tempera¬ 
tures and short operation durations. The coupling of these systems 
with building needs to be optimized in order to improve their 
energy, economic and environmental relevance. 

The control strategy is made outside the model by giving the 
mass flow rate, the inlet water temperature and the heat loads as 
inputs for the model (Appendix C). The mass flow can be fixed at 
a constant value or be variable to match a thermal output of the 
auxiliary, the water inlet temperature can be given by a water logic 
based on outdoor temperature and the heat loads can be given by a 
data file or a dynamic thermal simulation of a building. Its also pos¬ 
sible to implement “power led” or “cost led” control strategies by 
triggering the pCHP unit only if there is a sufficient power demand 
or if it is economically profitable: the heat can be stored in a buffer 
tank if there is no needs, otherwise, the heat demand can be covered 
by the auxiliary. 


7. Conclusion 

A semi-physical model of a pCHP device and a gas boiler has 
been developed and implemented in the TRNSYS environment in 
order to assess the energy, environmental and economic impacts 
of |jiCHP. The main advantages of the model are its simplicity while 
keeping a good accuracy. The model is able to simulate transient 
phases which can be important to assess the global efficiency of 
the device. 

The parameters have been identified by experimental investi¬ 
gations on a gas Stirling engine p,CHP unit. The initial method has 
been simplified in order to ease the identification of the model only 
with manufacturer data and a few experimental tests. The simpli¬ 
fied method does not degrade the accuracy of the model. The model 
has been validated for the Stirling engine for inlet cooling water 
temperatures from 30 to 60 °C and cooling water flow rates from 8 
to 161 min -1 . 

At the moment, this model is only adapted for gas applications: 
natural gas, propane or biogas. Its adaptation for biomass-CHP 
requires some improvements such as a CO-emission model during 
start-up/cooling phases or a new combustion model for the biomass 
fuel as Bouvenot et al. [34] did for a wood pellet pCHP device. 
Moreover, the model has been only validated for a non-modulating 
Stirling pCHP. Future works will be undertaken to extend the appli¬ 
cations of this pCHP model. 

Then, these models will be used in order to optimize their use in 
buildings by coupling them with buildings in numerical platforms 
like TRNSYS. The aim will be to identify the most suitable environ¬ 
ment and operation of these systems to maximize their efficiencies 
and their profitability. 
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Appendix A. Model parameters 

Tables A1-A3. 

Appendix B. Transient phases modelling approach and 
steady state to transient states transition conditions 

Figs. B1-B3. 
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Fig. Bl. Modelling approach of the start-up phase. 
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Table A1 

Model parameters of the Stirling engine. 



Steady state parameters 




Transient state parameters 


Parameter 

Value 

Unit 

Parameter 

Value 

Unit 

1 

mnom 

i,l cv 

16 

kgmn- 1 

21 

v start 

400 

s 

2 

pnom 

30 

°C 

22 

x p 

75 

s 


cw,i 




start 



3 

pnom 

fuel 

6640 

Wth 

23 

A C 

45 

s 

4 

our 

5315 

w th 

24 

A C 

75 

s 

5 

pnom 

1 gross 

864 

W e l 

25 

Atop 

250 

s 

6 

Pelectro 

4 

W e l 

26 


750 

s 

7 

Paux 

40 

Wel 

27 

A C 

85 

s 

8 

a 

12.35 

WthK- 1 

28 


12,350 

JK- 1 

9 

b 

-27.44 

W th (kgmn -1 ) _1 





10 

c 

-1.43 

w th K- 1 





11 

d 

8.62 

W th (kgmn- 1 )- 1 





12 

e 

-2.16 

W d K 1 





13 

f 

0.73 

W e i (kgmn- 1 )- 1 





14 

pnom 

exh 

38.8 

°C 





15 

g 

0.91 

- 





16 

h 

-0.31 

0 C(kgmn _1 ) _1 





17 

rh nom . 

cond 

2.16 E- 4 

kgs- 1 





18 

i 

2.84 E“ 6 

J kg- 1 °C“ 1 





19 

j 

-8.26 E- 6 

J kg- 1 (kgmn- 1 )- 1 





20 

[ UA ] i oss 

2.76 

WK- 1 






Table A2 

Model parameters of the additional boiler. 


Steady state parameters 




Transient state parameters 


Parameter 

Value 

Unit 

Parameter 

Value 

Unit 

1 

YY^nom 

IIL cv 

16 

kgmn -1 

18 

start aux 

450 

s 

2 

pnom 

30 

°C 

19 

At 0 - 

45 

s 


cv,i 




start aux 



3 

pnom 
fuel aux 

25,769 

w th 

20 

T° 

stop aux 

100 

s 

4 

r ) nom 

K AiXaux 

23,703 

Wth 

21 

At 0 

stop aux 

750 

s 

5 

Pelectro 

4 

Wei 

22 

[MCW 

12,350 

JK- 1 

6 

Paux 

60 

Wei 





7 

a' 

-256.1 

W th K- 1 





8 

b ’ 

594.14 

W th (kgmn- 1 )- 1 





9 

c’ 

-236.2 

w th K- 1 





10 

d' 

629.88 

W th (kgmn- 1 )- 1 





11 

j^nom 

38.2 

°C 





12 

g’ 

0.68 

- 





13 

h' 

0.08 

°C(kgmn“ 1 ) _1 





14 

fh nom 

cond aux 

7.58 E“ 4 

kgs- 1 





15 

r 

1.63 E- 5 

J kg- 1 °C“ 1 





16 

j ' 

-3.37 E“ 5 

J kg -1 (kgmn- 1 )- 1 





17 

[ IdA ] i oss 

2.76 

WK- 1 






Table A3 

Simplified model parameters of the Stirling engine. 


Steady state parameters 




Transient state parameters 


Parameter 

Value 

Unit 

Parameter 

Value 

Unit 

1 

nom 

,n cw 

16 

kgmn -1 

11 

‘'start 

400 

s 

2 

j'nom 

cw,i 

30 

°C 

12 

^start 

75 

s 

3 

pnom 

fuel 

6,749 

w th 

13 

A C 

45 

s 

4 

r\nom 
y HtX 

5,260 

w th 

14 

A ^start 

75 

s 

5 

pnom 

1 gross 

861 

Wel 

15 

Atop 

250 

s 

6 

Pelectro 

4 

Wel 

16 

A C 

750 

s 

7 

Paux 

40 

Wel 

17 

A C 

85 

s 

8 

a 

12.32 

WthK- 1 





9 

e 

-2.16 

W d K" 1 





10 

AT 

8 

°c 
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Fig. B2. Modelling approach of the cooling phase. 
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Fig. B3. steady state to transient states transition conditions. 


Appendix C. Algorithm chart 

Fig. Cl. 

- internal parameters surrounded in dotted orange line 

- 28 external parameters to be identified surrounded in solid 
orange line and written in black 


-17 external parameters to be identified of the simpli¬ 
fied model surrounded in solid orange line and written in 
violet 

- parameters which can be internal (default values) or external (to 
be identified) surrounded in solid red line 

- main outputs surrounded in blue line 

- main inputs surrounded in black line 
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Fig. Cl. Algorithm flow chart. 
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